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ABSTRACT 


The random number of failure-related events in certain repairable ageing systems, like certain nuclear power 
plant components, during a given time interval, may be often modelled by a compound Poisson distribution. 
One of these is the Polya-Aeppli distribution. The derivation of a stationary Polya-Aeppli distribution as a 
limiting distribution of rare events for stationary Bernouilli trials with first order Markov dependence is 
considered. But if the parameters of the Polya-Aeppli distribution are suitable time functions, we could expect 
that the resulting distribution would allow us to consider the distribution of failure-related events in an ageing 
system. If a critical number of damages produce an emergent failure, the abovementioned results can be applied 
in a reliability analysis. It is natural to ask under what conditions a Polya-Aeppli distribution could be a 
limiting distribution for non-homogeneous Bernouilli trials with first order Markov dependence. In this paper 
this problem is analyzed and possible applications of the obtained results to ageing or deteriorating nuclear 
power plant components are considered. The two traditional ways of modelling repairable systems in reliability 
theory: the “as bad as old” concept, that assumes that the replaced component is exactly under the same 
conditions as was the aged component before failure, and the “as good as new” concept, that assumes that the 
new component is under the same conditions of the replaced component when it was new, are briefly discussed 
in relation with the findings of the present work. 


1. INTRODUCTION 


The ultimate (but always elusive) goal of any nuclear engineer engaged in design, 
maintenance or operation activities related with nuclear power plants (NPP) would be to 
achieve a downtime as near zero as possible. A careful assessment of vulnerable areas of 
plant operation to allow the prevention of sequences of damage events that could result in the 
failure of some significant component is usually successfully done. Nevertheless, in some 
cases this is not feasible. A reliability approach must be applied to determine both the backup 
hardware needed in stock, and the average frequency and duration of repair operations. In 
these cases, quality assurance in design, production and operation allows the elimination of 
systematic or repeatable patterns of failure. An indication of achievement of a suitable 
optimization in these processes is the random nature of the events that cause failures and that 
the failures themselves are rare events. So, in this case the theory of stochastic processes is 
the appropriate framework for mathematical modelling. Within this framework the failure 
processes may be considered as chains of rare events or equivalently as point processes. For 


each stochastic process in continuous time with a discrete state space, there is a point process 
where the points are the random times of entry to a pre-assigned state or sets of states. In 
reliability theory what matters is the entrance in a failure state or the exit from it. 


The purpose of the present paper is to discuss some reliability aspects that could be of interest 
in the design and maintenance of NPP. The random number of damage or failure-related 
events M(t) during the interval (0,1) in certain repairable ageing systems may be composed, 


from time to time, by clumps of almost simultaneous rare events. If they are considered as 
truly simultaneous, it is possible to expect to be able to describe them by compound Poisson 
distributions (also known as “contagious” distributions) [1] defined in (0,r). It is conceivable 


that this same random number could be described also in discrete time, using random chains 
of rare events (failure- related) that are not independent from a probabilistic standpoint: we 
consider random damages or failures in such conditions that the occurrence of one of them 
modifies the probability of occurrence of the next one. Under suitable conditions, it could be 
expected that the probability distribution of the number of occurrences N, of a certain rare 


event in a great number of trials in a chain approaches a probability distribution of the 
compound Poisson type in continuous time, if the time interval Arbetween trials and the 
number of trials verify n-Ar=rwhile nis made to increase. The study of this kind of 


convergence of discrete distributions of failure-related events defined on chains, towards 
discrete distributions of failure-related events defined in continuous time, allows the 
interpretation of compound Poisson distributions and their related stochastic processes by 
models of more elementary nature in which well defined and meaningful mechanisms of 
random interaction may be introduced. In fact, Ar may be related with cyclical processes in the 
subsystems of the NPP. In certain reliability problems Arcan relate to the periods of 
mechanical vibrations induced in solid components of pumps, pipes, heat exchangers, fuel 
assemblies and structural parts of the nuclear isle. In other situations, it could have a 
connection with the slower periods of pressure and temperature oscillations. Thus, the 
parameters of the distributions of failure-related events may relate to the basic conditional 
probabilities of the chains. The following discussion is aimed to give a concrete example of 
the abovementioned method. It will be focussed in the Polya-Aeppli distribution of 
probabilities (also known as the Poisson-Geometric distribution), its relationship with random 
processes (both in the stationary and non-stationary cases), and its relations with Markov 
chains of rare events. 


If the repair time is negligible when compared with the time scale of the failure process, two 
approaches are often used in reliability theory in NPP: the non-stationary Poisson Process and 
the Renewal Process. As stressed by Saldanha, Frutuoso e Melo, and Noriega [2], the first 
approach may be too pessimistic, since it assumes that the system after repair is as bad as it 
was before repair and the second approach too optimistic, since it assumes that the system, 
after repair is as good as new. To get a balance between these traditional approaches, a more 
realistic modelling of failure and repair processes in an ageing system must be done. As 
noticed in [2], the explicit use of the complete intensity function allows a more realistic 
modelling of failure and repair processes in an ageing system, between the extremes of the 
“as bad as old” and the “as good as new” concepts. The theoretical approach proposed in the 
present work could be considered as a complement to the mathematical framework 
commonly applied in NPP reliability analysis. It is intended as a suggestion to further explore 
the feasibility of applying the approach in suitable NPP reliability problems. 
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2. POLYA-AEPPLI DISTRIBUTIONS, MARKOV CHAINS AND RELIABILITY 


2.1. The Polya-Aeppli distribution as limit of the distribution of rare events in a 
succession of Bernouilli trials with first order Markov dependence 


We first consider a succession of independent Bernouilli trials. Let p be the probability of 
Ein a single trial, andg the probability of E° (the complementary event, p+q=1) also ina 
single trial, and N,, the random number of events Fin n consecutive trials. If the mean value 
n-p Of N,,tends to a constant 4 whenn — 0, while the number 7 of trails increases without 
bounds, then it is well known that the probability P,(k)= P(N, =k), (k =0,1,2,...7) tends to a 


Poisson distribution. This is possible if p = p(n)= A Gite Y, with <«, >0 when n increases, so 
n 


that E becomes increasingly rare. The probability generating function (pgf) 


g,(z)= S"p, (k)-z* of the random variable N, is g,,(z)=(¢+ p-z)" = fi + (ten), (z ) 
k=0 n 


When nincreases, g,(z)approaches g(z)=e“""), the pgf of a Poisson distribution. A fairly 
general theorem of probability theory [3] allows us to conclude that if a sequence of 
pef g,(z)converges to a pgfg(z), the corresponding sequence P,(k)converges to P(k), the 
distribution of probabilities of the random variable characterized by g(z). Then, in fact, 
P,(k)converges to a Poisson distribution. 

Now, let us consider a succession of Bernouilli trials with first order Markov dependence (a 
Markov chain with two elements in its state space). In this case the probability of occurrence 
of E in the trial number n depends of the result of the trial number n—-1. Either directly as was 
done by Cernuschi and Saleme [4], or indirectly as was done by Bernstein [5], it is possible to 


k=n 
calculate the probability generating function (pgf) ¢,(z)= ¥,P,(k)-z* corresponding to the 
k=0 
random number of occurrences N,, of £ innconsecutive trials along the chain. The following 


formula is obtained ([4], [5]) g,(z)= (b(<)=Ag(z))-(” (<))" - (o(z)-4, (<))-(Ar(z))" (1) 


Ay(z)—Ag(z) 
The following equations are verified: 


a (z)= c(z)+ V(c(z))” -4-a(z) (2a) A,(z)= e(z)- V(c(z))? —4-a(z) (2b) 


2 2 
a(z)=(p'-p)-z (2c) b(z)=1-p-(I-z) (2d) e(z)=1-p+p'z  (2e) 
Here p is the probability of the event Fin a trial, if in the previous trial E“ occurred, and 


p'is the probability of Fif in the previous trial F occurred. 

The mean number of occurrences of Ein ntrials, w, is obtained from formula (1) for the pef, 

and the relation yz, = g', (1) (g’, (1) is the derivative of g,(z) taken forz=1). 

Thus we obtain: pea 2 p-(p-p) ( (p' »)') (3) 
1-(p'-p) 1+ pp! 

We can increase n without bound while maintaining “,, = 2 constant, if the probabilities 


ee 1 ae ; 
pand p' verify: p=pii)=—-(u+2,) (4a) p= P(n)= po +2,’ (4b) 
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Both ¢, and «', tends to zero while n tends tooo. Substituting (4a) and (4b) in formulae (2a) 

to (2e), from formula (1) it follows, after several fairly lengthy calculations ([6]), that g,, (z) 
(2-1) 

is asymptotically equivalent to (4,(z))” and both approach to: g(z)=e" Pe (6) 


(-pe)}z 
Equation (6) may be restated as a compound Poisson distribution: g(z)=e lees } (7) 
The distribution that corresponds to (7) is known as the Poisson-Geometric or Polya-Aeppli 
distribution. (If pg =0, (7) reduces to the pgf of Poisson distribution). Applying the same 
theorem used above in the case of independent Bernouilli trials, we conclude that P,(k), for 


Bernouilli trials of an increasingly rare event E with strong first order Markov dependence 
converges to a Polya-Aeppli distribution. 


2.2. Ageing systems, failure related rare events, and critical numbers to failure 
If the parameter w = 2-1, being 2a constant, and if pg is constant, then the distribution of the 
random number M(t) of failure-related events in (0,r) corresponds to a stationary process: a 


non ageing system. But if ~=[A(u)-dufor a certain A(r) (possibly increasing with time) and 
0 


Pe = Pa(t) (not necessarily increasing with time) the Polya-Aeppli distribution is non- 
stationary and could be used to describe an ageing system. 


If the critical number of failure-related events above which a failure is produced during (0,r) 
k=n, 

isn,, the reliability R(t) would be given by: R(t)= SY P(N(t)=4) (8) 
k=0 


It is possible to give a closed form analytical formulae for the probabilities P(N = k) = P(k) 
for the Polya-Aeppli distribution, but it is somewhat cumbersome to calculate directly with 
them ([7], [8]). If 1,is not too small (usually greater than 3 or 4) it is better to apply the 
method of steepest descent to estimate the probabilities P(k) fork >n,, as was done by 
Cernuschi and Castagnetto in their generalization of the distributions of Greenwood-Yule and 
Liiders [9]. In the case of the Polya-Aeppli distribution we obtain a fairly good 


approximation, beginning with the complex variable formula P(k)= 1 ate) -dz (where the 


oni zktl 


integral is taken on any simple closed contour that encloses the origin of the complex plane 


and excludes the singularity z= — of the function g(z) given by equations(6) or (7)). The 
P@ 


result of this procedure may be found in reference [10]. 


Then the reliability may be estimated from: R(t) = 1- FP (t)=k) 


k=n,+1 


2.3. Some remarks about non-stationary Polya-Aeppli distributions, non-homogeneous 
Markov chains and reliability in ageing systems 

It is natural to ask about a connection between non-homogeneous chains and non-stationary 
Polya-Aeppli distributions in the framework of reliability problems for ageing systems. 
Since the research of R. von Mises it is known that a sequence of independent Bernouilli 
trials with varying probabilities from trial to trial, gives a distribution for N,, that approaches 
to a Poisson distribution of parameter yz if the sum of the series of the probabilities of 
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jen 
occurrence of the event £ (that is the mean number of occurrences) E[N,,]= ¥ p(n) 
j=l 


converges to w and if max,<;<, (p j )tends to zero when the number 2 of trials tends to [11]. 


So asymptotically E is a rare event in every trial. A particularly interesting case, that may be 
generalized to Markov chains, is obtained if we assume a causative factor c(n)>1 such 


that p;,,(n)=c(n)- p(n) for every j. Now E[N,,]= Koy = - p,(n)so if this expression must tend 


to a constant, c(n)must tend to 1 from above and p,(n)must tend to zero as ms 
n 


The simplest model of a non-homogeneous sequence of Bernouilli trials with first order 
Markov dependence was constructed using the causative matrix concept by Harary, Lipstein 
and Styan [12]. Let us suppose that p',, = P(E,, /E,,;)is the probability of the E in the mtrial, 
given that E occurred in the m-1 trial, and q’,, = PE,’ i " Pm = PlE,, 1 E‘m-t 1; 

a,= P(E Smt Ema )the other conditional probabilities. With these conditional probabilities the 
probability transition matrix M,,,is constructed (1m, = P's. = Ps, ='m> Mon =F» )-Then 
if the chain has a constant causative matrix C, M 


mii =CM,, for every m=1,2,3,... So, if c, 
represent the elements of the causative matrix, for every m the following scalar equations 
apply: P' mst =C11* P'm*€12 °F m mst = 21° P’m*€22 °4'm3 

Pmt = C11 * Pm +12 ° Im Pmt = C11 * Pm +12 ° Im (9) 
The causative matrix is pseudo-stochastic in the sense that the off diagonal elements may be 
negative, but always c,,; +c). =1, c 3, +e =1 
In dealing with non-stationary chains it is advisable to work with vector probability 
generating functions that describe the distribution of the random number N,, of occurrences 
of E ina finite number of trials. So, let us introduce the (column) vector pgf g,,(z) with 
elements g',, (z)and g",, (z). The partial pgf g',, (z) describes the distribution of probabilities of 
Ein mtrials given that E occurs in trial number m, while g'',, (z)does the same thing when 
E does not occur in trial number m. 
Then: 8 mut (Z)= 2° m8 m (Z)+Z"Pm Bm (Z) 8" mat (Z)=4'n8'm (Z)+49m Bm’ (z) (10) 
In matrix terms: Z,,,,(z)=™,,,(z)g,,(z) From (9) and (10) this matrix equation is obtained: 
M nu (z)=C(z)M,,(z) . Being C(z) a matrix whose first row is the first row of C multiplied 
by z and whose second row is identical to the second row of C . From these equations the pgf 
for n consecutive trials g, (z)may be obtained in terms of powers of C(z)and M,. 
The spectral properties of the matrices M,,,(z)=C”"(z)M,(z)were already studied for the 


special case z =1by Le Maire and Gauffrey [13]. A possible extension of they results could 
be applied to discuss the asymptotic behavior of g,,(z)when bothn — » and the mean value of 


N,, converges. However, for our purpose here we give a heuristic argument closely related 


with the derivation of Polya-Aeppli distribution developed in 2.1.So, let us consider n trials in 
a non-stationary chain with transition probabilities from the trial number m to the trial m+1 
given by p'’,, (1), 4'm (2), Pm()s Im (n) (1<m<n) and the elements of the causative matrix are 


also functions of 1: c;;(2), ¢42(n),€2;(n), er (n) 

If E is increasingly rare asn >, so that q,,(n) > 1 (and as a consequence p,,(n) > 0), but the 
conditioning related to £ is so strong that p',, (n) > ps (and as a consequence 7’, (n) > 1- pe), 
then from equations (9) it follows that the causative matrix C(n) tends to the unit matrix. By a 
continuity argument we could expect that also for the case of a non-homogeneous sequence 
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of Bernouilli trials with first order Markov dependence and a causative matrix connecting the 
matrices of transition probabilities, the distribution of probabilities of the random number of 
occurrences of E inn trials converges to a Polya-Aeppli distribution. 


3. CONCLUSIONS 


(1) A connection between Markov chains and compound Poisson distributions, with focus in 
continuous random processes related with the Polya-Aeppli distribution that might be of 
interest from the standpoint of the mathematical modeling of some reliability problems posed 
by random mechanical vibrations and other oscillatory processes in NPP, was studied. 

(2) The approach developed in this paper may be considered of the conservative “as bad as 
old” type if the parameters of the asymptotic non-stationary Polya-Aeppli distribution 
describes an ageing system due to the cumulative effect of progressive damage processes. 

(3) Further research could be done to complete some of these findings and to extend the 
suggested approach to construct less-conservative models from the reliability standpoint. 
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